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Kinematic diffraction is well suited for a mathematical approach via measures, which has substantially been developed since 
the discovery of quasicrystals. The need for further insight emerged from the question of which distributions of matter, beyond 
perfect crystals, lead to pure point diffraction, hence to sharp Bragg peaks only. More recently, it has become apparent that one 
also has to study continuous diffraction in more detail, with a careful analysis of the different types of diffuse scattering involved. 
In this review, we summarise some key results, with particular emphasis on non-periodic structures. We choose an exposition on 
the basis of characteristic examples, while we refer to the existing literature for proofs and further details. 

1 Introduction 

Diffraction techniques have dominated the structure analysis _ ' . - . ° ' 

of solids for the last century, ever since von Laue and Bragg 
employed X-ray diffraction to determine the atomic struc- 
ture of crystalline materials. Despite the availability of di- 
rect imaging techniques such as electron and atomic force mi- 
croscopy, diffraction by X-rays, electrons and neutrons con- 
tinues to be the method of choice to detect order in the atomic 
arrangements of a substance; see Cowley's book 35 and refer- 
ences therein for background. 

In its full generality, the diffraction of a beam of X-rays, 
electrons or neutrons from a macroscopic piece of solid is a 
complicated physical process. It is the presence of inelas- 
tic and multiple scattering, prevalent particularly in electron 
diffraction, which makes it essentially impossible to arrive at 
a complete mathematical description of the process. Here, we 
restrict to kinematic diffraction in the far-field or Fraunhofer 
Jimit. In this case, powerful tools of harmonic analysis are 
available to attack the direct problem of calculating the (kine- 
matic) diffraction pattern of a given structure. 

In contrast, the inverse problem of determining a struc- 
ture from its diffraction intensities is extremely involved. A 
diffraction pattern rarely determines a structure uniquely, as 
there can be homometric structures sharing the same auto- 
correlation (and hence the same diffraction). 9 i 54 i 57 i 103 We are 
far away from a complete understanding of the homometry 
classes of structures, in particular if the diffraction spectrum 
contains continuous components. At present, a picture is 
emerging, based on the analysis of explicit examples, which 
highlight how large the homometry classes may be. 

Originally, much of the effort concentrated on the pure 
point part of diffraction, also called the Bragg diffraction, for 



Fig. 1: Experimental diffraction pattern of a quasicrystalline 
AlPdMn alloy. Figure courtesy of Conradin Beeli. 
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the case of ordinary (periodic) crystals, and later also for in- 
commensurate phases. Following the discovery of quasicrys- 
^ 70,79,92,1 18 w ifa tjjgjj- beautiful diffraction patterns, such as 
the one shown in Figure Q] a new mathematical approach was 
required. The associated paradigm shift also re-opened the 
discussion of what possible manifestations of order and dis- 
order in solids there are, and how these can be detected and 
quantified. While diffraction is one measure of order, the exis- 
tence of homometric structures of varying entropy 11-14 shows 
its limitations, as there are completely deterministic systems 
which cannot be distinguished from a randomly disordered 
system on the basis of pair correlations alone. Increasingly, 
the continuous or diffuse part of the diffraction is attracting 



attention, 



40.133.135 



not the least because improved experimen- 
tal techniques make the diffuse part accessible. Improving our 
understanding of diffuse diffraction is desirable, in particular 



1 




Fig. 2: A patch of the rhombic Penrose tiling. The arrow decorations 
of the edges encode the local rules. 

in view of the implications on disorder. 

The most successful approach to describe the structure of 
incommensurate crystals and quasicrystals employs additional 
dimensions. By embedding the ideal structure into a higher- 
dimensional 'superspace', it is possible to recover periodicity 
in the higher-dimensional space, and this picture can be ex- 
tended to cover certain aspects of random tilings as well. The 
standard tilings used to model the structure of quasicrystals 
are obtained in this way; for instance, the Penrose tiling 104 
shown in Figure|2]can be described as a projection 19 of a slice 
through the four-dimensional root lattice A4. Such structures, 
or their equivalent point sets, are called cut and project sets 
or model sets, and we shall discuss further examples below. 
Note that the Penrose tiling also possesses aperiodic, perfect 
local rules (or matching rules 52,75 ' 121 ), as well as an inflation 
symmetry. The local rules can be implemented as arrow dec- 
orations on the edges of the two rhombic prototiles, which, 
within any admissible patch, have to agree on all edges. These 
local rules are aperiodic in the sense that they are incompat- 
ible with any periodic tiling. They are perfect because they 
specify precisely the class of the rhombic Penrose tilings, in 
the sense that all space-filling tilings obeying these rules are 
locally indistinguishable (LI) from the rhombic Penrose tiling, 
the latter defined as a fixed point tiling of an inflation rule. 

It is worth noting that, while the lattice of periods of a peri- 
odic crystal is unique (though the choice of unit cell is not), 
there is considerable freedom in the choice of the building 
blocks of aperiodic tilings. In the case of the Penrose tiling, 



1 INTRODUCTION 




Fig. 3: The three allowed (pairwise) overlaps of the decagonal clus- 
ter. Overlapping markings are highlighted by colour. 

there exist a number of equivalent versions (in the sense 4 i 16 i 25 
of mutual local derivability), such as the Penrose pentagon 
tiling or the kite and dart tiling. One can even go beyond 
tilings and consider coverings of space. 80 In the case of the 
Penrose tiling, Gummelt's decagon covering 58 with a single 
cluster (and overlap rules encoded by the shading) has proved 
very popular, because it allows the description of a quasicrys- 
tal structure in terms of a single fundamental building block. 
The three allowed (pairwise) overlaps of the marked decagons, 
shown in Figure|3] are characterised by matching decorations. 
Figure |4] shows a patch of a corresponding covering, which is 
mutually locally derivable (MLD) with the Penrose tiling of 
Figure|2] 58,59 This covering also has an interpretation in terms 
of 'maxing rules', 49 ! 61 ' 72 where maximisation of one type of 
specified cluster leads to the Penrose rhombus tiling (up to 
zero density deviations). 72 Covering rules of either type have 
become quite fashionable in materials science. 109 i 124 For more 
examples on tilings, in particular on substitution tilings, we 
refer to the online Tilings Encyclopedia. 60 For the early de- 
velopment of the field, the reprint volume by Steinhardt and 
Ostlund^ 3 . is still a valuable source. 

This review attempts to present an overview of the develop- 
ment of mathematical diffraction theory in the 30 years since 
the discovery of quasicrystals by Shechtman et al.^& While 
we aim to provide the reader with a flavour of the mathemati- 
cal methods and assumptions, we will not dive deeply into the 
technical details. In particular, we will not present any for- 
mal proofs, though we do state several non-trivial results ex- 
plicitly. We refer to our recent review 13 and our forthcoming 
book, 16 and the references contained therein, for more details 
on the rigorous mathematical treatment. Three complemen- 
tary review volumes 3 ' 22 ' 97 with mathematical articles are also 
highly recommended. Here, we select examples that are both 
characteristic and somewhat supplementary to previous pre- 
sentations. 

In Section 12 we start with a concise summary of the sys- 
tematic approach using measures (in the mathematical sense, 
such as Lebesgue measure X, which is used to measure vol- 
ume in Euclidean space), which was pioneered in this context 
by Hof.— We first apply this approach to the diffraction 
of perfect crystals in Section [3] and then discuss the case of 
mathematical quasicrystals based on a cut and project scheme 
in Section |4] Like perfect (or idealised) crystals, these sys- 
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2 METHODS AND GENERAL RES ULTS 



2.1 Measures, convolutions and Fourier transforms 




Fig. 4: A patch of Gummelt's decagon covering. 

tems are pure point diffractive, which means that the diffrac- 
tion pattern consists of sharp (Bragg) peaks only. Afterwards, 
in Section [5] we proceed to systems with continuous diffrac- 
tion, covering both the case of singular continuous and abso- 
lutely continuous diffraction by means of representative exam- 
ples, including a probabilistic model for thermal fluctuations. 
In particular, we consider random tilings, which are relevant 
because most quasicrystalline materials show entropic stabil- 
isation and therefore are expected to include configuration al 
disorder. 

2 Methods and general results 

For a satisfying mathematical approach, we should exclude 
any boundary effects, and hence consider infinite systems that 
represent the scattering medium. Traditionally, there are two 
seemingly contradictory ways to describe a system, either in 
terms of functions which represent the density of the scat- 
tering medium, or by lattices or, more generally, tilings of 
space, whose decorations mimic the atomic positions. This 
dichotomy has sparked some rather fierce disputes between 
the tiling school and the density function school, in particular 
in the years following the discovery of quasicrystals. How- 
ever, the two viewpoints can be reconciled by embedding them 
into a more general frame. One way of doing that is to intro- 
duce measures, which comprise (almost) periodic functions 
and tilings as special cases. As measures quantify distribu- 



tions in spaces, this approach is in fact very natural, and well 
suited to describe both the distribution of matter in the scat- 
tering medium and the distribution of (scattered) intensity in 
space. We therefore start by briefly introducing the concepts 
and main properties that will be needed in our context. 

2.1 Measures, convolutions and Fourier transforms 

Due to the Riesz-Markov representation theorem,i2& it is pos- 
sible to think of a measure as a linear functional, i.e., as a 
linear map that associates a number to each function from an 
appropriate space. A (complex) measure on W 1 is then a lin- 
ear functional (with values in the complex numbers C) on the 
space C c (K d ) of complex- valued, continuous (test) functions 
of compact support, subject to the condition that, for every 
compact set K C M d , there is a constant a K such that 

\n(g)\ < flrlUH- 

for all test functions g with support in K. Here, WgW^ = 
sup vGA -|g(x)| is the supremum norm of g. 

We write fl(g) or J M d g(x)dji(x) for the measure of a func- 
tion g, and jJ.(A) — /^(1 A ) for the measure of a set A C M. d , 
where 

Jl, ifxeA, 
1 0, otherwise, 

denotes the characteristic function of the set A. 

If jl is a complex measure, the conjugate of ii is the measure 
fl which is defined by g H> /J.(g). A measure is called real (or 
signed), when jl = jl, and it is called positive when ji(g) > 
for all g > 0. For every measure fl, there is a smallest positive 
measure, denoted by \fi\, such that |jU(g)| < for all non- 

negative g. This is called the total variation (or absolute value) 
of fl. A measure jl is called finite or bounded, if = 
|jU | (W 1 ) is finite, otherwise it is called unbounded. As we want 
to describe infinite point sets in space, we usually deal with the 
latter case, but we will assume that measures are translation 
bounded. This means that, for any compact set K c Mr, the 
total variation satisfies 

sup \n\(t+K) < e», 

fGR rf 

so wherever you move your compact set K, its total variation 
measure is always finite. 

2.2 Autocorrelation and diffraction measures 

If A C Mr is a point set that is a Delone set (a set where points 
neither get arbitrarily close nor so sparse that it accommodates 
arbitrarily large empty balls), the corresponding Dirac comb 34 

8a := 

.VGA 
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3 DIFFRACTION OF PERFECT CRYSTALS 



is a translation bounded measure, where 8 X is the normalised 
(Dirac) point measure at x (so 8 x (g) — g(x), or, in the for- 
mal notation used in physics, J R( / g(y) 8 (y — x) dy — g(x)). In 
what follows, we use such Dirac combs to represent the scat- 
tering medium, possibly with (in general complex) scattering 
weights w(x) at position x g M. d . The corresponding weighted 
Dirac comb is denoted as 



ft) 



.VGA 



If ft) is a translation bounded measure, the corresponding 
diffraction measure is the Fourier transform of the autocorre- 
lation measure, where we shall assume that the latter exists. 
In any given example, this has to be verified, of course. The 
autocorrelation measure of ft) is defined as the limit 



7 = 7 ffl = ft) © ft) 



:— lim 



a\R*a\R 

VOl(fifl) 



(1) 



where Br denotes the open ball of radius R around € M.' 1 . By 
co\ R we denote the restriction of ft) to the ball Br. For a mea- 
sure pi, its 'flipped-over' version pi is defined via pi(g) = ji(g), 
where g(x) ~ g{—x). The operation * is the ordinary convo- 
lution of measures, which is a generalisation of the standard 
convolution of integrable functions, 



(f*g)(x) := [ f(x-y)g(y)dy= [ f(y)g(x- 
For finite measures pt and v on Mr, it is defined by 



y)dy. 



(m*v)(s) 



g(x + y)dn(x)dv(y) 



for any function g g C C (R ), which is then again a finite mea- 
sure. The volume-averaged convolution © (also called the 
Eberlein convolution, in analogy to a similar approach— in 
the theory of almost periodic measures) is needed in Eq. (Q3, 
because ft) itself is generally an unbounded measure and the 
direct convolution is not defined. For example, if A denotes 
the standard Lebesgue measure (for volume), A * A is not de- 
fined, while A © A = A . 

If the autocorrelation measure 7 of ft) exists, its Fourier 
transform 7 does as well, and 7 is a translation bounded, pos- 
itive measure, called the diffraction measure of ft). It corre- 
sponds to the kinematic scattering intensity observed in an 
experiment in the sense that it quantifies how much scatter- 
ing intensity reaches a given volume in li-space. Relative to 
Lebesgue measure A, the diffraction measure has a unique de- 
compositionJiS 

y Ipp ' Ysc ' Tac 

into its pure point part (the Bragg peaks, of which there are 
at most countably many), its absolutely continuous part (the 
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Fig. 5: Illustration of the distribution function F(x) of the classic 
middle-thirds Cantor set. The iterative construction for the latter is 
sketched in the inset. 



diffuse background scattering, which has a locally integrable 
density relative to A) and its singular continuous part (which 
simply means anything that remains, which is nothing in many 
standard cases considered in crystallography). Each of the 
three terms is again a positive measure. Singular continu- 
ous measures are weird objects: they give no weight to single 
points, but are still concentrated to an (uncountable!) set of 
zero Lebesgue measure. A well-known example is the proba- 
bility measure for the classic middle-thirds Cantor set,^S& with 
the 'Devil's stair case' as its distribution function, which is 
constant almost everywhere; see Figure |5] Singular continu- 
ous diffraction does occur in realistic models though, — and 
should not be disregarded. 

3 Diffraction of perfect crystals 

In our setting, a perfect (infinite) crystal in af-space is a lattice- 
periodic (discrete) structure. It is defined by its lattice of pe- 
riods T C W l and the decoration of a fundamental domain of 
T, which together completely specify the distribution of scat- 
terers in space. It is therefore described by a crystallographic 
measure 

ft) = fi*8 r , (2) 

where ji is a finite measure. The latter can be chosen as the 
restriction of ft) to a fundamental domain of F. Depending on 
the nature of fi, the resulting measure ft) can be pure point or 
continuous (for instance, if jj. is the constant measure on the 
fundamental domain, ft) would be proportional to Lebesgue 
measure), or a mixture of both types. One can think of the 
Dirac comb 8 r as implementing the lattice periodicity, while 
pi describes the distribution of scatterers in a fundamental do- 
main of r. 

The autocorrelation of the crystallographic measure ft) of 
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3.2 Diffraction of crystallographic structures 



Eq. © is given by 

7= dens(r) (jU *ju) * 8 r , (3) 

which follows by using the relation 8 r = 8 r together with 
8 r © 8 r = dens(F) 8 r . Here, dens(F) denotes the density (per 
unit volume) of the lattice T, which is the reciprocal of the 
volume of its fundamental domain. Consequently, y is also 
a T -periodic measure. In order to obtain the corresponding 
diffraction measure, we need to know how to calculate the 
Fourier transform of lattice-periodic measures. 

3.1 Poisson's summation formula 

A powerful tool for the Fourier analysis of lattice-periodic 
measures is the Poisson summation formula (PSF). For a lat- 
tice r cM. d (which means that T is a discrete subgroup of 
M. d such that the factor group M. /T is compact), the Fourier 
transform of the corresponding Dirac comb 8r is 

8^. = dens(r) S n , (4) 

where F* denotes the dual or reciprocal lattice of T. The 
latter is defined by 

T* = {x € R d I (x\y) £ Z for all y € F} . 

Here and below, (x\y) denotes the scalar product of x,y S M. d . 
Note that sometimes a factor 2n is included in the definition 
of the reciprocal lattice, which we prefer to incorporate in our 
definition of the Fourier transform. For a suitable function <$>, 
our convention for Fourier transform is 

£(jfc) / e- 271 '^ (p{x)dx, 

where k,x € W 1 and again (k\x) denotes their scalar prod- 
uct. The Fourier transform y of a positive definite measure 
y (which means that Y(g*g) > holds for all g e C c (E. d )) 
is defined as the unique extension 29 108 of the Fourier trans- 
form of functions. It is conveniently defined in the setting of 
tempered distributions, 108 which provide concrete means to 
calculate the transforms. 

By the Bochner-Schwartz theorem, 108 the diffraction mea- 
sure is then a translation bounded positive measure. In addi- 
tion, we will make use of the convolution theorem for mea- 
sures. This states that if jj. is a finite measure and v a transla- 
tion bounded measure on R rf , the convolution fx * v exists and 
is a translation bounded measure. 29 If V is not only a tempered 
distribution, but itself also a measure, one has the convolution 
identity ji * v =fiv. The latter is then again a measure, which 
is absolutely continuous relative to v, because /i is a bounded, 
uniformly continuous function on K. in this case. 



3.2 Diffraction of crystallographic structures 

Using the PSF together with the convolution theorem, the 
Fourier transform of the crystallographic autocorrelatio n mea- 
sure y of Eq. (f3]l can be calculated as 

y=(dens(r)) 2 |£| 2 5 r *. (5) 

Clearly, this is a pure point measure, concentrated on the 
dual lattice T*. Note that ki is a uniformly continuous and 
bounded function that is evaluated only at points of the dual 
lattice T* . While different admissible choices for the measure 
il (describing the same system) lead to different such func- 
tions, they agree on all points of F*, so that the result does 
not depend on this choice. If y{{k}) = for some k € F*, one 
calls this an extinction. Extinctions are characteristic features 
of further symmetries, also of generalised type. 

3.3 Planar a-phases 

Let us consider an interesting example in some detail. Starting 
from a checker board, viewed as a decoration of the square 
lattice, we assume that the grey squares are stiff (or solid), 
while the white squares are empty. One can now twist the 
structure by rotating the grey squares alternately in opposite 
directions by an angle <p € (— f , f ), 




This way, a new periodic structure emerges where the white 
squares are deformed into congruent rhombuses. This struc- 
ture is the lattice-periodic repetition of the motif above, and 
resembles a planar ff-phase and related quasicrystal approx- 
imants.— A couple of examples are shown in Figure [6] The 
second is related to structures found in 12-fold symmetric qua- 
sicrystals.— 

We consider the associated Dirac comb 

<*V = 8 R v s * S a<pZ 2 • 

obtained by placing a normalised point (or Dirac) measure 
at each vertex point. Here, we have = 2cos(<p) and 

R <P = (Z(l) "oS?)' While S = {°' e i' e 2> e i + e i\ denotes 
the vertex set of the unit square [0, l] 2 . The corresponding 
diffraction measure is obtained via Eq. §5^ as 

1 +cos(27r(^ (p e 1 |fc)) 1 +cos(27t(R (p e 2 \k)) 

Y f = 2cOs(<p) 2 2cOs(<p) 2 5 Z2/2cos(,p)' 



5 



4 DIFFRACTION OF MATHEMATICAL QUASICRYSTALS 




• ■ • ■ ■ • • • • • ■ • 

• •< ••• • ■ • • • . 

• • • • • • • 

-•• ••• 

• • • • - • • 

. . • . * . ••• ■ . • 

• ••• • • •■■ 

• •••••• • 

• • • • • • • • 

• •• . • * . - • • • ••• 

• •••••• • 

• ••• •■• ••• 

. . • ■ m ■ ••• • • • 




■ • • 
• • • 



• ■ • 



• • • 

• • • 



• ■ • 

• • • 



• • • ■ 

• • • • 

• • • • 

• • • ■ 



Fig. 6: Planar rj-phases with angles <p = n/S (top) and <p = k/12 
(bottom), shown with the correct relative length scale. In the latter 
case, the rhombus dissects into two equilateral triangles. 



with (x\y) denoting the scalar product in M 2 . 

When <p = (which means we are back to the square lat- 
tice), this expression reduces to f = 8 z2 , as it must, while 
inserting <p = ±7r/4 leads to f^Z = 45 fl which reflects 

jr/4 

the double weight of the point measures at each vertex in this 
limit. For angles <p with tan(<p) irrational, one has extinctions 
precisely for all wave vectors k — ('p-,'—) with m,m 1 = and 



Fig. 7: Diffraction patterns for the two rj-phases of Figure [6] All 
distances and intensities are shown in the correct relative scale. 



tion and that has an area proportional to the intensity. This 
choice resembles the experimental situation in a reasonable 
way. Both patterns are non-periodic, due to the incommensu- 
rate positions of the points in the fundamental cell. While all 
Bragg peaks are located at positions of the corresponding dual 
lattices, there is an apparent approximate 8- or 12-fold sym- 
metry in the patterns (sometimes called pseudo-symmetry), 
which is why we chose these examples. 



e 2Z + 1 . When tan(p) is rational, there are further 4 Diffraction of mathematical quasicrystals 



extinctions, which can be calculated from the explicit formula 
for the diffraction measure y^. 

The diffraction patterns for the two examples (<p = tt/8 and 
<p = n/\2) from Figure [6] are illustrated in Figure A Bragg 
peak is represented by a dot that is centred at the peak posi- 



We now leave the realm of lattice periodic systems to discuss 
aperiodically ordered structures, in particular quasicrystals. 
Before we move on to structures with non-crystallographic 
symmetries, let us briefly consider the inclusion of incommen- 
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4.1 Incommensurate phases 



surability in a lattice periodic system, which can be seen as a 
first step towards the structure of mathematical quasicrystals. 

4.1 Incommensurate phases 

The systematic investigation of incommensurate systems was 
pioneered by de Wolff 42 and by Janner and Janssen. 71 We re- 
fer to a recent monograph by van Smaaleni 2 ^ and references 
contained therein for details and background, and concentrate 
on a couple of elementary examples here. 

The simplest incommensurate structure arises from com- 
bining two periodic Dirac combs with incommensurate peri- 
ods, such as 

<»« ■= 8 z + 8 a . 

with a > irrational. While this is unphysical in the sense 
that positions of scatterers become arbitrarily close, it is in- 
structive to look at the diffraction for this toy model. Observe 
the Eberlein convolutions 8 Z © 8 aZ = ^A, which is a conse- 
quence of a being irrational, and 8 aZ © 8 aZ — jjS aZ , which 
follows from a simple density calculation. Then, the autocor- 
relation turns out to be 

1 2 
Ya = 8 Z + - S aZ + -^, 

which leads to the diffraction measure 

Ya = °Z+tf 8 Z/a + - 8 a 

by an application of the PSF together with A = 5 . This 
pure point diffraction measure reflects the two periodic con- 
stituents. There are Bragg peaks on the integer lattice (with 
intensity 1) and on the reciprocal lattice Z/a of the lattice 
ali, with intensity a~ 2 . Note that the intensity of the central 
peak is 1 + a~ 2 + 2a~' = (1 +a~') 2 , in line with the density 
of the underlying point set. One might expect that the relative 
position of the two constituent lattices does not matter, which 
indeed is the case. Introducing a relative shift u between the 
two periodic combs does not affect the result, in the sense that 
the diffraction of the Dirac comb co a u = 8 Z + 8 ll+aZ is still 
given by y a , independently of the value of u. 

While this system is of limited practical relevance in one di- 
mension, one can build higher-dimensional systems using the 
same idea. This results in incommensurate systems which are 
called composite structures. Let us discuss a simple example. 
Fix some a > and consider the Dirac comb 

(0 = 8 r _ + 8„ . r = 8 z2 + 8„ * 8, , 

where F = aZ x Z C K 2 is a planar lattice, and u <E 1R 2 an 
arbitrary shift. For a <G Q, the underlying point set is crystal- 
lographic, with Z 2 nF as its lattice of periods. Here, we are 



oooooooo 
oooooooo 
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oooooooo 
oooooooo 
oooooooo 



Fig. 8: Composite structure comprising atoms on the square lattice 
(black dots) and on the shifted lattice u + T (circles), with shift u = 
\\,\) and lattice r = aZxZfora = T= ^(l + \/J). 

interested in the non-periodic case, so let us assume that a is 
irrational. An example is displayed in Figure [8] 

The autocorrelation for the Dirac comb CO evaluates as 

Y = 5 Z 2 + ^8 r + 1 (8 U + 8_ U )*(X <g> 8 Z ) , 

where fl ® V stands for the (tensor) product of two measures. 
The Fourier transform of y can be obtained by applying the 
Poisson summation formula and the convolution theorem. It 
has the form 

1 2 

Y = 8 i? + ^2 8 r* + - cos(2^ 2 m 2 ) (S ® 8 Z ) 

with the dual (reciprocal) lattice T* = (^Z) x Z. Note that 
the final term only involves the second components of k and 
m, due to the presence of the term 8 in the measure (so only 
k i = contributes). In the diffraction measure, the composite 
structure is visible via additional intensities of the peaks along 
the vertical axis. The total intensity of a Bragg peak at position 
(0,«) with n € Zis 

7 ({(0,»)}) = 1 + -L + 1 cos(27T«« 2 ) > (1 - I) 2 > 0. 

The corresponding diffraction pattern for the example of Fig- 
ure |8]is shown in Figure [9] 

Of course, this is merely a sketch of any real system. For a 
more realistic system, one should take into account the mod- 
ulation in the positions induced by the different local neigh- 
bourhoods. 119 . 128 . 131 . 132 
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4 DIFFRACTION OF MATHEMATICAL QUASICRYSTALS 



Fig. 9: Diffraction pattern of the composite structure of Figure [8] 
Each Bragg peak is again represented by a dot which is centred at the 
position of the peak and whose area is proportional to the intensity. 
One can clearly recognise the peaks on the two lattices Z and F* , and 
the alternating intensity of the peaks along the vertical axis, which are 
due to the choice u 2 = j ■ 



Here, we consider a simpler case, based on the modulation 
of a periodic structure. A modulated structure arises by locally 
displacing positions of a crystalline point set, ensuring a min- 
imal distance between the new positions. For example, start 
with the integer lattice Z and deform it by moving the points 
according to a real-valued displacement function h. The de- 
formed point set is then given by 



A h = {n + h(n) | n € Z}, 



(6) 



and 8 Af denotes the corresponding Dirac comb. To be 
concrete, consider the displacement function h(n) = e{an}, 
where a and £ are real numbers and where {x} = x — [x] de- 
notes the fractional part of x. Since \h(n)\ < £, the deformed 
point set respects a minimum distance between points, as long 
as £ is sufficiently small. Clearly, if a is a rational number, 
the resulting point set is once again periodic, while it is non- 
periodic for irrational values of a, which is the case we are 
interested in here. 

To understand the corresponding set A h , it is advantageous 
to use an embedding in the plane, known as the 'superspace 
approach' in crystallography. 128 Define a planar lattice as the 
integer span of two basis vectors 



1 

-a 



where we use the notation (u,v) z = {mu+nv \ m, n £ Z}. Con- 
sider now the line pattern obtained as the F-orbit of the line 




Fig. 10: Superspace approach for the modulated point set A h of 
Eq. for e = 0.35 and a « 0.2941. The lines (or 'targets') inter- 
secting the horizontal axis are shown in black. 



from the origin to the point (£,1) (with the end point not in- 
cluded). Then, A h is the set of intersections of the horizontal 
axis with these line segments; see Figure [K)]for an illustration. 

Using the fact that, for irrational a, the sequence of num- 
bers ({ocn}) neZ is uniformly distributed in the unit interval, 85 
one can calculate the autocorrelation y h of the Dirac comb on 
A h explicitly. The result is 



Yi, 



meZ 



{am}) 8, 



m+e{cem} 



+ {am}8 m _ £(l _ 



e(\-{am})) 



The corresponding diffraction measure y, reads 

% = £ \A(k)\ 2 8 kl 

keZ[a] 

with (complex) amplitudes 

A(k) = e- m * sinc(jik*) , 



(7) 



(8) 



where sinc(x) = sin(x) fx. The map k i-> k* acts on elements of 
Z[a] = {r + sa \ r,s e Z} as (r + sa) ^ (re + s(l + ea)) for 
any r,s 6 Z. In this example, % is a pure point measure which 
is supported on a dense set. Despite the denseness of the Bragg 
peaks, the total intensity scattered into any compact subset of 
R is finite, because the intensities are locally summable. The 
proof for the diffraction formula is non-trivial. However, this 
can be interpreted as a special case of the diffraction of model 
sets (cut and project sets), because the modulated structure © 
is in fact a model set. We now turn our attention to this general 
notion, and discuss a number of relevant examples and their 
diffraction. 
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4.2 Model sets 

There are a number of ways to construct aperiodically ordered 
systems. 55 From the viewpoint of diffraction, the best under- 
stood is a natural generalisation of lattice-periodic structures 
obtained by a projection from a higher-dimensional lattice. 
Such systems are called cut and project sets or model sets, 98 
and can be produced in a number of essentially equivalent 
ways, 50 including de Bruijn's grid method— and Kramer's 
'Klotz construction',— as well as a number of other ap- 
proaches. 66 ' 99 

The model set approach can be viewed as a generalisation 
of the notion of a quasiperiodic function.— In the simplest 
setting, the idea is much like what we saw for the modulated 
phase in Figure [10] above: The aperiodic structure emerges 
by taking a cut across a higher-dimensional periodic structure, 
using a direction that is incommensurate with the lattice. The 
general setting for the case of Euclidean model sets is encoded 
in the cut and project scheme (CPS) 



l-i 



U 



U dense 



(9) 



model atomic structures of quasicrystals. However, as long 
as there exists a local rule to switch from the point set to the 
tiling picture and vice versa, we can consider both structures 
as equivalent (as any atomic structure will be a local decora- 
tion of either), or shortly as MLD (which stands for mutual 
local derivability). 4 For instance, in one dimension, a tiling of 
K by two intervals of different lengths is clearly MLD with the 
set of left endpoints of all intervals. 

In what follows, we only consider regular model sets, so 
we require that the boundary dW of the window W has zero 
Lebesgue measure. The Euclidean setting (O generalises to 
the case where the internal space is a locally compact Abelian 
group. 96 i 98 ' 115 We shall meet an example later, where the in- 
ternal space is based on 2-adic integers. 

Regular model sets are pure point diffractive, 23 i 6 V 15 and 
in this sense are natural generalisations of lattices. This is 
a central result of the theory of model sets which has been 
proved by methods of dynamical systems theory, 67 i 9 M 15 in 
terms of almost periodic measures 23 ' 101 ' 126 and, following a 
suggestion by Lagarias, by using the Poisson summation for- 
mula for the embedding lattice and Weyl's lemma on uniform 
distribution. 16 The diffraction measure y of the Dirac comb 
8 A is explicitly given by 



where W 1 is the physical (sometimes also called direct or par- 
allel) space, and W" is referred to as the internal (or perpendic- 
ular) space. Here, «£? c W !+m is a lattice in d + m dimensions, 
and n and n in( denote the natural projections onto the phys- 
ical and internal spaces. It is assumed that L = 7t(-S?) C W 1 
is a bijective image of ££ in direct space, and that the set 
L* = 7T int (Jz?) C M!" is dense in internal space. As a conse- 
quence, the *-map 98 x h- > x* is well-defined on L. 
A model set for a given CPS is then a set of the form 



A = {xeL\x* ew} : 



(10) 



where W C K m (called the window or acceptance domain) 
is a relatively compact subset of W" with non-empty inte- 
rior. More generally, also translates of such sets are called 
model sets. The elements of the model set A lie in the pro- 
jected lattice L in direct space, and the window in internal 
space determines which elements of L are selected. The con- 
ditions on the window ensure that the model set A is a Delone 
set. In fact, a model set A is always a Meyer set, 96 98 which 
means that A — A := {x — y | x,y £ A} is uniformly discrete, 
while A is relatively dense. Note that uniform discreteness 
of A — A implies that of A , and is actually a much stronger 

condition. 86,87,9^98 

Clearly, the projection approach produces point sets in 
space rather than the tilings that are conventionally used to 



: E \A(k)\ 2 

kei® 



(11) 



Here, L® = 7f(_S?*) is the corresponding Fourier module, 
which is the projection of the higher-dimensional dual lattice. 
The amplitudes are given be 



_ dens(A) p 
A{k > ~ vol(W) lw 



(-**), 



(12) 



where l w is the characteristic function of the window W. Var- 
ious generalisations, in particular to certain weighted Dirac 
combs, have been discussed in the literature, 16 i 23 - m - 115 An al- 
ternative (and somewhat complementary) approach based on 
an average periodic structure can be employed to unravel var- 
ious modulation features in the diffraction patterns of qua- 
sicrystals. This is systematically explained in a recent re- 



,136 

for further details. 



\ ie\\ by Wolny and coworkers; see references cited there 



4.3 One-dimensional examples 

We start by re-expressing the modulated point set A/, of 
Eq. ( [Tol l as a cut and project set. To this end, we need to write 
Af, via an orthogonal projection, rather than via the (implicit) 
skew projection of Figure[l0] This can be done by introducing 
the matrix A = (q j £ ) and considering the lattice Jz? = AT. 
This lattice and its dual lattice are given in terms of a Z-basis 
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Fig. 11: Model set description of the modulated point set A/, of 
Figure [Tol 



by 

if: 



-a 



-e 
1 



if* 



a 

1 +ea 



The two generating vectors and the lattice points of if are 
shown in Figure [TT] 

The set A/, is now a model set for the CPS with lattice 



if C R = 



(so d = m = 1 and both direct and inter- 



nal space are R). The window is the interval W = [0, — 1), and 
the condition x* £ W selects all lattice points that are located 
within the shaded strip of Figure QT](which is the reason why 
this approach is sometimes also referred to as the strip projec- 
tion method). For e = 0, we get a (non-minimal) embedding 
of Z in R 2 , and for rational a = | with coprime integers p and 
q we obtain a periodic point set with lattice of periods qL. 

The formulas (|7]i for the diffraction and (|8) for the ampli- 
tudes now follow from the general result of Eqs. ([TT1 and (flZt . 
The Fourier module is L® = 71: (if*) —Z[a], and the action of 
the *-map can be read off from the explicit bases of if and 
if* given above. 

The most frequently invoked example of a one-dimensional 
(mathematical) quasicrystal is the Fibonacci chain. Its geo- 
metric version is built from two intervals (prototiles) L and 
S (for long and short) of lengths x = (1 + V5)/2 and 1. It 
can be generated by iterating the square of the inflation rule 
L 1 — ^ LS, S M> L, starting from a legal seed (such as L\L, where 
the vertical line indicates the reference point). This leads to 
the two-sided interval sequence 

... LSLLSLSLLSLLSLSLLSLSL \ LSLLSLSLLSLLSLSLLSLSL ... 

The bi-infinite sequence is aperiodic, with relative frequencies 
and T~ 2 for the two prototiles. 




Fig. 12: Model set description of the Fibonacci chain. 

Define two point sets Al and As as the left endpoints of the 
corresponding intervals in the chain, taking the reference point 
as 0. They are model sets for the CPS © with d = m = 1 and 

L = Z[t] = {m + nx | m,n G Z} . 

The corresponding planar lattice is 

3f = 



1 — T 



which has density 1 / V 5 and the dual lattice 
2t-1 /Yt-1\ ( 1 



J^* 



One has Al $ — {x E L\x* £ W^s} with the windows 



W L = (-l,T-2] and 

and the *-map defined by V5 n- 
m + n — nx. The construction is illustrated in Figure [12] The 
Fibonacci model set is A = Al U A5, with window 



W s = (t-2,t-1] 
-V5, so that (m + m) 



W = W L UW S = (-1,T-1]. 

This way, A is a point set of density t/ \/5 = (t + 2)/5. Note 
that it is possible to modify the embedding lattice if by scal- 
ing the internal space relative to physical space. In particular, 
multiplying the scale of internal space by t, the embedding 
lattice is a rotated copy of y/t + 2 1?. 

The Dirac comb co = 8j\ is pure point diffractive, by an 
application of the model set diffraction theore m 23 ' 67 ! 115 men- 
tioned before. The diffraction measure 7 is explicitly given by 
Eq. (fTTb with the amplitudes 



A(k) 



(13) 



via Eq. i\2i , where sinc(x) = sin(x)/x. The phase factor re- 
flects the position of the window, which is centred at (t — 2)/2. 
The sum in Eq. ( fTTb runs over the Fourier module 

L® = 7T(if*) = ^=Z[T]. 
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5 10 15 20 

Fig. 13: Diffraction pattern for the Fibonacci chain A . The Bragg 
peak at has height (dens(A)) 2 = (t+1)/Sw 0.5206, and the entire 
pattern is reflection symmetric. 

A sketch of the diffraction pattern is shown in Figure [13] Note 
that the intensity function I(k) = |A(£)| 2 vanishes if and only 
if Tk* € Z \ {0}. This means k = ix with O^eZ. Since all 
such points lie in the Fourier module L®, we have identified 
all extinctions. These are a fingerprint of the intrinsic inflation 
symmetry. 

As an example of a limit-periodic structure, consider the 
period doubling sequence. Written as an element w G {0, 1} Z , 
it is given by w(2n) = 1, w(4n + 1) = and w(4n +3) = w(n) 
for n G Z. This rule specifies every position except n = — 1, 
where we can choose either possibility. Both possibilities can 
also be obtained as a fixed point sequence of the square of the 
substitution 1 n- 10, 11. The two sequences have cores 

...101 11010101 1101^1 101 11010101 11011... 

and are locally indistinguishable. They thus define the same 

system. The underlying Toeplitz structure of a hierarchy 

of scaled and shifted copies of Z is apparent from the for- 
mulai^l 

A = {neZ\w(n) = l} = (J ((2-4^Z + (4 £ - 1)) 

e>o 

for w(— 1 ) = (with — 1 added to A in the other case). This set 
can be described as a model set, but with the internal space be- 
ing the 2-adic integers. Consequently, the diffraction measure 
of the Dirac comb 5^ is again pure point. 

The corresponding diffraction formula can be given explic- 
itly as follows. 12,16 The Fourier module is 

L® = Z[|] = {f | (r = 0,meZ) or (r > l,modd)}, 

so that we can again use Eq. ( fTTT i. Here, the amplitudes are 





.,,.,.1.,,.,. 










( 


) 









Fig. 14: Absolute values of the diffraction amplitudes for the period 
doubling chain. The diffraction pattern is 1-periodic. 



with i:= | el*. This parametrisation specifies k uniquely. 
Figure [14] shows the absolute values \A(k)\ for k G L® n [0, 1]. 
This pattern repeats Z-periodically. 

Further one-dimensional examples will be discussed in Sec- 
tion |5] in the context of continuous diffraction measures. Let 
us now turn our attention to higher-dimensional model sets. 



4.4 Cyclotomic model sets 

For the description of two-dimensional tilings, it is advanta- 
geous to work with complex numbers x + iy in C rather than 
with points (x,y) in M 2 . In C, a rotation by an angle <p just 
corresponds to multiplication with the complex number e"? ' . 
This point of view is a natural generalisation of de Bruijn's 
method 41 and the Fourier space approach.— A natural way 
to implement an n-fold rotational symmetry is to choose a 
primitive nth root of unity 4, G C (so 4" = 1 an d 4™ 7^ 1 f° r 
1 <m<n), and to consider the Z-module Z[<i;J of cyclotomic 
integers, comprising all integer linear combinations of powers 
of ^„ (the solutions of the equation x" = 1). One can think of 
cyclotomic integers as the set of all points in the plane that can 
be reached by taking steps of unit length along the directions 
of a regular n-star. Clearly, the resulting point set is symmetric 
under rotations of multiples of 2n/n; in fact, under rotations 
by multiples of n/nifn is odd. Therefore, one usually restricts 
to integers n ^ 2 mod 4 to avoid duplications. 

The cases n G {1,2} are trivial in the sense that the resulting 
point sets lie on the real axis. The choices n G {3,4} lead to 
crystallographic point sets, the triangular lattice with sixfold 
symmetry and the square lattice with fourfold symmetry. Any 
other choice n > 5, n ^ 2 mod 4, produces a dense point set 
in the plane, with «-fold symmetry for even n, and 2n-fold 
symmetry for odd n. 

The dense point set Z[<i;J can be embedded into a lattice by 
lifting it to a suitable higher-dimensional space, essentially by 
making all directions in the n-star that are linearly indepen- 
dent over the integers (there are (n) such directions, where 
is Euler's totient function) also linearly independent over the 
real numbers. A natural way to do this is the Minkowski (or 
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Fig. 15: Ammann-Beenker tiling as a cyclotomic model set. 



Galois) embedding 

_£f„ = {(x,G 2 (x 



which defines a lattice Jz?„ C 



dx))\xeZ[Q} (14) 



Here, o>, with 



I < £ < 0(h), are the Galois automorphisms of the corre- 
sponding cyclotomic number field, mapping a primitive root 
% n H> %n ' to a primitive root % n ' , where {mi | 1 < i < 0(n)} = 
{1 < A; < n | and n coprime}, together with a suitable order- 
ing. Note that is the identity map. Using the lattice J£„ in a 
cut and project scheme, with physical space R 2 ~ C and inter- 
nal space M^"' -2 , we produce cyclotomic model sets, which, 
for suitably chosen windows, have n-fold (2«-fold) rotational 
symmetry. 

As an explicit example, we consider the classic Ammann- 
Beenker (or octagonal) tiling 2 ' 28 as a cyclotomic model set 
with n = 8. Other standard examples of this type include 
the Penrose tilingi^ 4 - of Figure [2] and the Tubingen trian- 
gle tiling 19 (both with tenfold symmetry) and Gahler's shield 
tilin g 47 ' 48 (with twelvefold symmetry). The latter is locally 
equivalent (MLD) with a tiling introduced by Socolar.i2£ 
Since 0(5) = 0(8) = 0(12) — 4, all these tilings are obtained 
from cut and project schemes (O with internal space R 2 . 

Of course, the resulting tilings are only rotationally sym- 
metric if the window is chosen to have an appropriate ro- 
tational symmetry. To obtain the (undecorated) Ammann- 
Beenker tiling, the window W AB has to be chosen as a regular 
octagon, of unit edge length. The module 

L = = {«0+«1^8+«2^8 +"3^8 I (no,ni,n 2 , ni ) eZ 4 } 

is dense in the plane, and naturally lifts to a hypercubic lat- 
tice in four dimensions (the corresponding Minkowski embed- 
ding Jdg is a scaled and rotated version of Z 4 ). The *-map 



■ • • 
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Fig. 16: Diffraction pattern of the Ammann-Beenker tiling. 



can be chosen as the Galois automorphism <^ 8 n- and the 
Ammann-Beenker model set is then obtained as 



A 



AB 



{xeL\x*eW AB }. 



Figure [15] shows the picture in physical and internal space. 
Selecting points x E L whose *-image falls inside the octag- 
onal window (shown on the right of Figure IT3T > produces the 
point set in physical space shown on the left. Connecting all 
points of unit distance (which clearly is a local rule) recovers 
the Ammann-Beenker tiling, which is MLD with the cyclo- 
tomic model set. The decorations needed for the approach via 
local rules 2 add some non-local information, and cannot be 
recovered from the undecorated tiling alone. 48 i 120 . 

The diffraction of the Dirac comb on the Ammann-Beenker 
model set can be calculated via Eqs. (fTTT > and (fT~2t . It is a pure 
point measure supported on the dual module L® = (with 
the factor j due to the aforementioned scaling of the hypercu- 
bic lattice in the Minkowski embedding). The amplitudes (or 
Fourier-Bohr coefficients) are 

because the lattice J?? 8 has density j . A central patch of the 
diffraction image, obtained via an exact calculation of the 
Fourier transform of the octagonal window, is shown in Fig- 
ure Q2] In principle, the diffraction of any model set can be 
calculated (at least approximately) in this way, although it may 
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be complicated if the window is not a simple polygon or cir- 
cle, such as for the square-triangle tilings where the windows 
have fractal boundaries. 18 i 63 

4.5 Icosahedral model sets 

The model set approach works in any dimension. In particu- 
lar, it can be used to construct icosahedrally symmetric tilings 
in three-dimensional space, which are particularly relevant for 
applications in crystallography. The minimum embedding di- 
mension for this purpose is six, because one needs a faithful 
action of the icosahedral group and an invariant subspace of 
dimension 3. In this setting, there exist three different classes 
of icosahedral model sets, which correspond to the three dif- 
ferent hypercubic lattices (primitive, face-centred and body- 
centred) in six dimensions. 112 i 116 As body-centred icosahe- 
dral structures have not yet been identified in quasicrystals, we 
concentrate on the other two classes, and discuss one example 
of either type. 

For the primitive icosahedral tiling, we start from a lattice 
Jz? that is similar to the integer lattice Z 6 , and use a cut and 
project scheme (0 where both physical and internal space are 
R 3 . The corresponding window is shown in Figure [17] it is a 
semi-regular polyhedron known as Kepler's triacontahedron. 
The triacontahedron has edge length V2 + T, volume 20t 3 and 
surface area60f, where T = (1 +V5) /2 is again the golden ra- 
tio. This approach was pioneered by Kramer and Neri,— and 
the tiling is also sometimes called the Ammann-Kramer-Neri 
tiling (Ammann described the tiling earlier by different means, 
without publishing his findings; compare the corresponding 
comments in Mackay's early paper 93 ). Some authors also call 
it the three-dimensional Penrose tiling, in analogy to the five- 
fold rhombus tiling in the plane. 

The primitive icosahedral tiling is built from two rhombo- 
hedral prototiles, a thick (or prolate, called T v ) and a thin (or 
oblate, called T ) one. They can be defined as the convex hulls 
of their vertices 

T p = conv{0 ) v 1 ,v 2 ,v 3 ,v 1 +v 2 ,v l + v 3 ,v 2 + v 3 ,v 1 + v 2 + v 3 }, 
T o = conv{0,v l ,v 2 ,v 5 ,v 1 +v 2 ,v 1 +v 5 ,v 2 + v 5 ,v l +v 2 + v 5 }, 

where the basis vectors are^^ 

V 1 = (T,0,1), V 2 =(T,0,-1), V 3 = (1,T,0), 
V 4 =(-1,T,0), V 5 = (0,1,T), V 6 = (0,-1,T). 

These six vectors generate the primitive icosahedral module 

= (vi,V2,V3,v 4 ,v 5 ,v 6 ) Z) 

which plays the role of L = 7C(Jf) in the corresponding cut and 
project scheme ©. The *-map acts as (a,b,c) >-> x{a' ,b',c') 
on ^# p , where ' denotes algebraic conjugation (which maps 




Fig. 17: Kepler's triacontahedron as window of the primitive icosa- 
hedral tiling due to Kramer and Neri.— 

Vj H> -Vj, hence %' = 1 — t). In this formulation, the em- 
bedding lattice Jz? = {(x,x*) \ x g L} is similar to Z 6 , and ex- 
plicitly generated by the Z-basis {(v ; ,v*) | 1 < i < 6} with 
the vectors from Eq. (fl5t . Consequently, the fundamental 
cell of Jz? has volume 40(4t + 3), so that the density of Jz? 
is (7-4t)/200. 

A sketch of the two prototiles is shown in Figure [18] The 
rhombohedra have solid angles n/5, 3n/5 and In/ 5 as in- 
dicated. The solid angles in both cases add up to An. The 
prototiles have volumes 2t 2 (for T p ) and 2t (for T ). Note 
that ten rhombohedra of each type can be assembled to 
fill Kepler's triacontahedron of Figure [17] 

Figure [19] shows the only vertex star out of the 24 possible 
vertex stars of the Kramer-Neri tiling which has full icosahe- 
dral symmetry. In any tiling obtained from a generic model 
set, this vertex type occupies a subset that itself is a model set 
with the T -3 -scaled triacontahedron as its window. This prop- 
erty corresponds to the invariance of the module under 
multiplication by T 3 and reflects the inflation symmetry of the 
primitive icosahedral tiling. The corresponding (local) infla- 
tion rule, however, turns out to be rather complicated and has 
never been presented in complete detail. 




Fig. 18: Sketch of the two rhombohedral prototiles. 
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Fig. 19: Icosahedrally symmetric vertex star of the Kramer-Neri 
tiling, comprising 20 acute rhombohedra. 
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The diffraction of the Dirac comb on the primitive icosahe- 
dral model set can be calculated again by Eqs. (fTTT i and ( fT2] >. 
The Fourier module in this case is 



f P ® = 



1 



2(t + 2) 



'■? ■ 



The diffraction spectrum consists of a dense set of Bragg 
peaks located on L®, of which only a discrete subset has inten- 
sity above any chosen (positive) threshold. A full calculation 
of the Fourier transform of the triacontahedron was given by 
Elser, 45 so the intensities can be obtained explicitly. 

For simplicity, however, we employ a spherical approxima- 
tion to the amplitudes, by replacing the triacontahedron by a 
sphere of equal volume 20t 3 . The radius of the sphere turns 
out to be 

/ 15 \ 1/3 
R = [— J T fa 2.7246. 

Because the triacontahedral window is well approximated by 
this sphere, the difference between the approximate and the 
exact diffraction intensities is tiny, and irrelevant for our pur- 
pose. Note that the approximation only affects the values of 
the amplitudes, not the location of the peaks (except for ex- 
tinctions, which might show up in the approximation as tiny 
intensities). The Fourier transform of the spherical window 
evaluates as 



vo\(B R 



i Jbr 



2Mk*y 



dy 



3(sin(z) -zcos(z)) 



with z — 27t\k*\R. Figure l20l shows sections through the cor- 
responding three-dimensional diffraction patterns, orthogonal 
to the fivefold, threefold and twofold symmetry axes. 

An example of an F-type (face-centred) icosahedral model 
set is Danzer's tiling^ which was first constructed from an 
inflation rule, and is also known as the ABCK tiling, after the 
labels Danzer used for the four tetrahedral prototiles. In the 



• . - •: 
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Fig. 20: Fivefold (top), threefold (middle) and twofold (bottom) 
sections of the diffraction pattern of the primitive icosahedral tiling. 
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Fig. 21: The tiles of the Danzer tiling appear in groups of four 
(A,B,C) or eight (K), forming (topological) octahedra. 




Fig. 22: The three icosahedrally symmetric vertex stars of the 
Danzer tiling, comprising exclusively tiles of type B, C or K. 

ABCK tiling, the tetrahedral tiles always occur in the configu- 
rations shown in Figure[2T] so one can alternatively work with 
assembled prototiles consisting of four tiles of type A, B or C, 
and eight tiles of type K. 

The ABCK tiling is mutually locally derivable 39 ' 113 from 
the Socolar-Steinhardt tiling,i2£ so both describe equivalent 
structures. An interesting property of Danzer's ABCK tiling is 
the fact that it possesses particularly simple perfect local rules, 
which can be formulated as purely geometric packing rules 
on the level of the octahedra. 38 The Danzer tiling has three 
icosahedrally symmetric vertex stars, each comprising just one 
type of tiles, which are shown in Figure [22] Under inflation, 
these act as seeds of globally icosahedrally symmetric Danzer 
tilings. 

For the Danzer tiling comprising these larger prototiles of 
Figure[2T| all vertices are located on the face-centred icosahe- 
dral module 

= <V 1 + V2,V2 + V3,V 3 -|-V4,V4+V 5 ,V5-r-V 6 ,V 6 -V 1 > z , 

which is a sub module of of index 2. Explicitly, one has 

= My U {y#f + ? 2 m) , 

where u = j(vj — v 2 + v 3 — v 4 + v 5 — v 6 ) = (1,1,1). For this 
choice of coordinates, u is not in The vertex point set can 
be described as a three-component model set 81 ' 113 based on 
a cut and project scheme (0 with physical and internal space 
R 3 . The corresponding lattice Jz? is the embedding of in 
R 6 , which is similar to the root lattice Ds- The vertices of the 
four types of (topological) octahedra (thus disregarding their 
centres) separate into three different types, which stem from 
different cosets of the embedding lattice. 




Fig. 23: Windows for the vertices of type I, II and III of the Danzer 
tiling. They are shown in the correct relative size and orientation. 

In fact, the usual description as a three-component model 
set uses the projections of so-called 'holes' in the lattice Jz? . 
Holes are vertices of the Voronoi cells whose distance from 
points of the lattice is a local maximum.— If the distance is an 
absolute maximum, the hole is called deep, otherwise shallow. 
The vertices of the Danzer tiling then fall into three groups: 
Vertices of type I are projections from deep holes which lie 
in the coset jSf + (1, 1, 1,T,T, t), those of type II from deep 
holes in the coset j£f + (t, t, T, — 1, — 1, — 1) and vertices of 
type III from shallow holes in the coset _S? + (t,0, 1, — 1, 0, t). 
The corresponding three windows have icosahedral symme- 
try and are shown in Figure [23] The window for vertex type 

1 is a dodecahedral extension of an icosahedron, with pen- 
tagonal edge length 2 and volume 20(4 — t), the window for 
vertex type II is a dodecahedron of edge length 2/t and vol- 
ume 4(t + 2), and the third window is a great dodecahedron 
(a Kepler-Poinsot polyhedron), with pentagonal edge length 

2 and volume 20(t— 1). The *-map is the same as for the 
primitive model set above. 

The diffraction pattern of the Danzer tiling has spots on the 
corresponding dual module 

^f® = (-^p U {JK V + u)) (16) 

with u as above. Whereas the primitive tilings has diffraction 
spots on = 2(x+2) ^P on ly> me Danzer tiling has addi- 
tional spots on the shifted copy 2 (x+2) (*^P M )' Note that the 
union „# p U (.<# P + u) = ./# B corresponds to the body-centred 
icosahedral module. 

Due to the relation between the symmetry directions and the 
shift u, not all high-symmetry sections through the origin will 
show peaks from both modules in Eq. (fl6l l. In fact, only the 
twofold sections through the origin contain peaks from both 
parts in Eq. ( fT6l and thus display the full Fourier module, 
while the three- and fivefold sections only contain peaks from 
2{z+2) ^P' FigureGUsfiows the twofold section of the diffrac- 
tion for a Dirac comb of vertex type II only, so the window 
is simply a dodecahedron, which we approximate by a sphere 
of radius R = (^^1) 1/3 « 1.5118. In Figure[24] the 'black' 
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Fig. 24: Twofold section of the diffraction pattern of the ABCK 
tiling, with scatterers on all vertices of type II. 

diffraction spots belong to 2R+2) ' / ^P' wri ^ e tne § re y s P ots be- 
long to the coset. 

To visualise the diffraction along the fivefold axis, we com- 
bine the section through the origin with two parallel sections 
containing the coset reflections. The result is shown in Fig- 
ure!^ The spots in 2 (t+2) ^p aie a 8 am shown in black, while 
the two different grey colours distinguish the spots in the two 
parallel planes containing u (dark grey) or — u (light grey). 
This pattern demonstrates that the overall rotational symmetry 
here is fivefold (not tenfold) and inversion symmetric. The lat- 
ter property accounts for the tenfold rotation symmetry of the 
section through the origin (black spots). Sections with three- 
fold symmetry display the analogous phenomena. 

The distinction between the diffraction of a primitive and 
of a face-centred icosahedral model set is thus immediately 
recognisable from the spot locations in a twofold section. For 
further (practical) details and examples we refer to the recent 
literature J2£ 

Within the realm of regular model sets, diffraction is thus 
pretty well understood. We know that regular model sets are 
pure point diffractive, 23 i 67 i 1 15 and Eqs. ( fTTT i and (Y2\ provide 
explicit expressions for the intensities in terms of the Fourier 
transform of the window. Homometry of model sets (within 
the same cut and project scheme) can be traced back to equal- 
ity of the covariogram of the window, and explicit examples 
of homometric model sets have been constructed. 9 Thermal 
fluctuations can be taken into account in a fashion that is anal- 
ogous to the crystallographic case; 6 - 68 see Section l531 below. 

The Bragg diffraction has some robustness property beyond 
the class of regular model sets. Recently, Strungam 126 proved 
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Fig. 25: Fivefold section of the diffraction pattern of the ABCK 
tiling; see text for details. 

that, for any Meyer set A C Mr, the corresponding Dirac comb 
CO = 8 A always shows a non-trivial point diffraction, though in 
general the spectrum will be mixed and not pure point. How- 
ever, the point part is substantial in the sense that for any 
e > 0, the set of peaks [k G W 1 | Y({k}) > (1 - e) 7({0})} 
(all peaks with intensity near the maximum intensity) is rela- 
tively dense. While we do not have a complete answer to the 
question what structures are pure point diffractive, 32 it is clear 
that a pure point spectrum imposes strong constraints on the 
possible structures.— 

For the remainder of this article, we are looking at systems 
that show continuous diffraction, both singular and absolutely 
continuous. The discussion of examples with and without ran- 
dom disorder will shed some light on the much more complex 
situation beyond the pure point diffractive regime. 

5 Systems with continuous diffraction 

It seems a relatively recent experimental observation that 
diffuse scattering (as an indication of structural disorder, 
and not just of thermal fluctuations) is a widespread phe- 
nomenon. 13 V 35 It is thus natural to also investigate continu- 
ous diffraction spectra from a more mathematical perspective. 
Again, we briefly present illustrative examples, most of which 
have been analysed completely and rigorously by now. 

5.1 Singular continuous diffraction 

Let us begin by recalling the paradigm of singular continuous 
diffraction, the Thue-Morse (or Pruhet-Thue-Morse) system. 1 
It is usually defined via the substitution rule a i-> ab, b i-> ba. A 
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Fig. 26: Distribution function of the Thue-Morse diffraction mea- 
sure on the unit interval. 



bi-infinite fixed point sequence w emerges from iterating the 
square of this rule with the legal seed a\a. Define the Dirac 
comb 

where f(a) — 1 and f(b) = — 1, One can now show that the 
autocorrelation measure exists 10 - 73 - 94 - 134 and is of the form 



7 = E iH^ 



with r\ (0) = 1 and the recursion 

r](2m) = Tj(m) and rj(2m+l) = — ^ (Tj(m) + Tj(m+1)), 

which is valid for all m € Z. This exact renormalisation-type 
structure is the golden key to prove the spectral type and to 
calculate the measure explicitly. 

The diffraction measure is 1 -periodic, 5,8 and hence of the 
form 7 = ji * 8 Z with a positive, singular continuous measure 
jX. To describe the latter explicitly, one defines the distribution 
function F(x) = fj,{[0,x]) on the unit interval. It is consistently 
extended to a function on R by setting F(x + n) = F(x) + n 
for n E Z. This way, F(x) — x is 1-periodic and possesses the 
uniformly converging Fourier series 

v-i I ?( OT ) 

Fix)— x= ) sm(27tmx) . 

w ^ mn ' 

m= 1 

For computational purposes, however, it is advantageous to 
use an approximation in terms of a uniformly converging se- 
quence of distribution functions as follows. Define F (x) = x 




Fig. 27: Patch of the squiral tiling, obtained by two inflation steps 
from the central seed, which is legal. 



and the functional iteration 

1 f 2x 

w*) = 2y < 



1 -cos(ny))dF N (y) 



for N > 0. Since this iteration maps distribution functions 
for absolutely continuous measures to distribution functions of 
the same type, one can write dF N (x) — f N (x) dx with a Radon- 
Nikodym density f N . One can now check explicitly that this 
leads to 



i(2'xx)), 



where the empty product is to be evaluated as 1. Since the 
densities f N become increasingly spiky (and do not converge 
as a sequence of functions), one uses the distribution functions 
F N to illustrate the resulting measure. Note that the sequence 
(F N ) Ne ^ converges uniformly, 8 but not absolutely. This is in 
line with the fact that fi is singular continuous, and thus can- 
not be approximated by a norm-converging sequence of abso- 
lutely continuous measures. 108 The resulting limit distribution 
function F is illustrated in Figure [26] Despite its similarity 
with the Cantor measure of Figure [5] F is a strictly increasing 
function. This means that there is no proper plateau here. 

A non-trivial planar example emerges from the squiral in- 
flation rule from Figure 10.1.4 in Griinbaum and Shephard. 56 
It effectively leads to an aperiodic 2-colouring of the square 
lattice, according to the chirality of the square dissections; see 
Figure 127] for an illustration. Positioning a point measure of 
weight 1 or — 1 in the centre of the two types of squares, one 
obtains a weighted Dirac comb with average weight 0. Due to 
the inflation structure, one can derive a recursion formula for 
the corresponding autocorrelation. 15 ' 16 

By constructive methods, in complete analogy to the case of 
the Thue-Morse sequence, one can show that this Dirac comb 
leads to a purely singular continuous diffraction measure. 15 ' 16 
As in the one-dimensional case, it can explicitly be calculated, 
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Fig. 28: Third step of the Riesz product structure of Eq. (1171 for the 
diffraction measure of the squiral tiling; see text for details. 



and represented as a two-dimensional Riesz product. The re 
suit reads 

N-l 

f N (*,y) = n*( 3 '*' 3 *>0> 

where the function # is defined by 



(17) 



&(x,y) = — yl +2cos(27Dc) + 2cos(2jiy) 

— 2cos(27r(x + y)) — 2cos(27r(x — y)) 



As in the one-dimensional case, the corresponding distribution 
function possesses a uniformly convergent Fourier series rep- 
resentation, which involves the autocorrelation coefficients. 
The density function / 3 (bottom) and the corresponding dis- 
tribution function F 3 (top, normalised such that F 3 (Q,Q) = 0) 
are shown in Figure [28] 



5.2 Rudin-Shapiro chain and its Bernoullisation 

A simple, deterministic system with continuous diffraction is 
the binary Rudin-Shapiro chain. While it is usually presented 
via a four-letter substitution rule, the corresponding weighted 
Dirac comb o RS = EneZ w ( ,I )5 n can be defined by the se- 
quence of weights (w(n))„ eZ with w(n) £ {±1}, initial condi- 
tions w(— 1) = — 1, w(0) = 1, and the recursion 



for £e {0,1} 
for £e {2,3} 

The arrangement of the two weights looks as follows 



»(4»+/)=r (B) ' — (is. 

1 ' \(-l)"+ l w(n), for £e {2,3}. 



• ■•0**0«f*tOOO«OOOOtOO«OM*OOOfoj«tOt»0««*«000«OM«0»«0*000»«fOt 

where the line denotes the origin, and filled (open) dots corre- 
spond to weights 1 (—1). 

Despite the deterministic structure, the autocorrelation mea- 
sure of the balanced Dirac comb co RS (which has average scat- 
tering strength 0) can be shown 105 i 114 i 117 to be y RS = 8 . A 
simple proof of this is obtained by considering the induced re- 
cursion relation for the autocorrelation coefficients. U ' 13 The 
corresponding diffraction measure is thus y^ s = X, which is 
purely absolutely continuous, and shows no trace whatsoever 
of the underlying deterministic order in the sequence. In fact, 
the system is (almost surely, meaning for almost all realisa- 
tions of the random sequence) homometric with the random 
Dirac comb on Z with weights from {±1} chosen at random, 
independently at each position and with equal probability. 

We can combine the deterministic sequence and indepen- 
dently chosen random numbers by considering the Bernoulli- 
sation of the Dirac comb a^g, which we define as 11 



a P = H w{n)X(n)8 n . 



(19) 



Here, (w(n)) n€Z is the binary Rudin-Shapiro sequence of 
weights from Eq. (fT8l l. whereas (X(n)) neZ is an i.i.d. family 
of random numbers, each taking values 1 and — 1 with prob- 
abilities p and 1 — p (so < p < 1), respectively. The limit 
cases p g {0,1} bring us back to the deterministic Dirac comb 
±0^5 , while the case p — i corresponds to the Bernoulli comb 
with weights 1 and —1 mentioned above. The Bernoullisation 
thus interpolates between the deterministic Rudin-Shapiro se- 
quence and the completely uncorrected sequence of indepen- 
dent random numbers. It can also be interpreted as a 'model 
of second thoughts', where the sign of the weight at position 
n is changed with probability 1 — p. 

Using the strong law of large numbers,— it can be 
shown 1116 that the autocorrelation y p of the Dirac comb co p 
is almost surely given by 



Y P = CP - I) 2 7rs + 4/7(1 - P) S = 5, 
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irrespective of the value of the parameter p £ [0,1]. So the 
diffraction of this Dirac comb, for any choice of the parameter 
p, is (almost surely) y p = X, and the entire family of Dirac 
combs is homometric. 

This simple example highlights the fact that diffraction in 
general cannot distinguish 'order' in the sense of a determinis- 
tic structure from that in the presence of entropy. Note that the 
deterministic Rudin-Shapiro sequence has zero entropy, while 
the Bernoulli comb has entropy log (2), which is the maximum 
entropy for a binary sequence. For general p, the entropy is 
H(p) = —plog(p) — (1 — p)log(l — p), so it varies continu- 
ously between and log(2). Regardless, the diffraction of all 
these combs is the same. This result provides a glimpse at 
how degenerate, and hence difficult, the inverse problem can 
be in the presence of continuous spectra. Similar arguments 
can be used in higher dimensions (in particular by considering 
product structures), and examples in two dimensions involv- 
ing lower rank entropy have also been discussed. 14 i 27 

5.3 Random displacements and thermal fluctuations 

There are various important applications of Bernoulli-type 
disorder in real systems. The most obvious one is known 
as the random occupation model, which covers lattice gases 
and models of chemical disorder. Traditionally, this has been 
formulated for lattice-based systems only, 35 ' 130 but the cor- 
responding results hold in much greater generality. This in- 
cludes model sets, 21 but also structures with a substantial de- 
gree of positional disorder. 6 i 83 i 84 It turns out that the lattice 
assumption can be replaced by rather general principles from 
probability theory that revolve around the strong law of large 
numbers. 46 

This change of perspective is also of value for the treat- 
ment of the effects of thermal fluctuations to the diffraction 
of solids. In fact, rather than restricting to a lattice and small 
vibrations in a harmonic potential, the famous Debye-Waller 
contribution 35 can alternatively be derived from the assump- 
tion that the scatterers are randomly displaced from their equi- 
librium positions, independently of each other, but based on 
the same probability distribution. This opens the door to an- 
other application of the strong law of large numbers, as was 
first observed by Hof. 68 Two further advantages are the valid- 
ity for considerably more general point sets than lattices and 
the independence of the argument of the small displacement 
assumption. At least for sufficiently high temperatures, this 
alternative approach is reasonable. 

Consider a Delone set A C Mr 1 that is sufficiently nice 
(where we refer to the literature 6 ' 68 for the precise conditions). 
In particular, we assume that the Dirac comb 8 A possesses the 
autocorrelation y. The random displacement is described as 

A' = {x + t x \x£A}, 



where {t x ) xeA is a family of i.i.d. random translation vectors 
with common probability distribution v. Then, with probabil- 
ity one, 8 A , has the autocorrelation 

y' = y* (v* v) + dens(A) (8 - v* v) . (20) 

The corresponding diffraction is obtained by Fourier trans- 
form and reads 

?' = |v| 2 y+dens(A)(l-|v| 2 ). (21) 

Here, visa uniformly continuous function on Mr that vanishes 
at infinity, and the formula holds almost surely, as Eq. (l20l . 
If y, the diffraction of 8 A , is a pure point measure, the pure 
point part of 7' is given by |v| 2 7 (hence by a modulation of 
the intensities, which is the Debye-Waller factor), while the 
continuous part is dens(A)(5 — V * v). Note, however, that 
Eq. (Bil l is by no means restricted to pure point diffractive 
systems. An explicit dependence on the temperature can be 
modelled by the appropriate choice of the displacement distri- 
bution v. Further details and generalisations are discussed in 
the literature .^ 3 - 

5.4 Random tilings 

Random tilings form a particularly interesting and relevant 
class of structures, as was early pointed out by Elser.— The 
structure of the various ensembles and their diffraction is not 
as well understood as in the deterministic case, though a fairly 
complete picture was sketched by Henley.- 6 " 2 From a physical 
point of view, most results are 'clear', on the basis of convinc- 
ing (scaling) arguments from statistical mechanics. The math- 
ematical counterpart, however, is still incomplete, and various 
properties have escaped a proof so far, particularly in dimen- 
sions 2 and higher. In fact, it is a characteristic feature of ran- 
dom tilings to show a strong dependence on the dimension, as 
we will illustrate by some examples. 

Let us first consider a random version of the Fibonacci 
chain. Here, one starts with two prototiles as before (one inter- 
val of length T and one of length 1 ), and builds a tiling of R by 
choosing them with probabilities p and 1 — p, where p = T~ 1 
leads (almost surely) to realisations with the same relative tile 
frequencies as the deterministic chain of Figure [12] Due to the 
linear arrangement, the ensemble is well under control by el- 
ementary methods from probability theory. In particular, one 
can either invoke the ergodicity of the Bernoulli (coin toss- 
ing) chain 17 or the renewal theorem 6 to show that the ran- 
dom Dirac comb obtained this way almost surely leads to the 
diffraction measure 

y = 8 Q + h(k)X (22) 
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Fig. 29: Continuous part of the diffraction pattern of a Fibonacci 
random tiling. The range for the wave number k on the horizontal 
axis is the same as in Figure [T3l 

with the Radon-Nikodym density function 

h(k) = — (sin(^/T)) 2 

5 T 2 (sin(7T/tT)) 2 + T(sin(7i;/t)) 2 -(sin(7rA:/T)) 2 ' 

The factor (t+2)/5 0.7236 is the density of the correspond- 
ing point set, which equals that of the deterministic counter- 
part discussed earlier. Apart from the trivial Bragg peak at 
k = 0, the diffraction is thus absolutely continuous. Figure [29] 
shows the function h, which is smooth but still shows a spiky 
structure that resembles the pure point diffraction of the per- 
fectly ordered Fibonacci chain from Figure Q~3] to an amazing 
degree. 

The mechanism behind the absolutely continuous nature of 
the diffraction in Eq. (l22l can be understood as follows. Due 
to the choice of the intervals, each realisation can be lifted 
within the cut and project scheme of the perfect Fibonacci 
chain of FigureQ~2] Almost surely, one then obtains a sequence 
of lattice points that deviate from the perfect case via fluc- 
tuations that diverge linearly with the system size. 62 ' 64 This 
destroys the coherence needed for Bragg peaks (at k ^ 0) or 
singular continuous contributions to y. 

Random tilings in the plane show a different behaviour, 
which also depends on the symmetry. In particular, it is impor- 
tant whether one deals with a crystallographic symmetry (such 
as statistical three- or sixfold symmetry in the lozenge tiling) 
or not (such as statistical eightfold symmetry in the random 
octagonal tiling). An example of the former case, with broken 
symmetry, is illustrated in Figure [30] The underlying ensem- 
ble is well studied in statistical physics. 43 i 64 i 74 i 76 i 77 

The lozenge (or rhombus) with opening angle n/3 occurs 
in three possible orientations in all typical lozenge random 
tilings (which are subject to the condition that any resulting 
tiling is face to face and covers the plane without overlaps). 
One can now use the relative frequencies of the three prototiles 
to parametrise the ensemble. By purely group theoretic meth- 




Fig. 30: Typical patch of a rhombus (or lozenge) random tiling, 
with periodic boundary conditions. Here, the vertical rhombus is less 
frequent than the other two types, hence breaking the statistical three- 
fold symmetry. 

ods, one can then show that the entropy has a unique maxi- 
mum at the (unique) point of maximal symmetry. 62,110 This 
shows an interesting entropic mechanism for the stabilisation 
of tilings with statistical symmetry. The value of the entropy 
(calculated per tile) is known exactly from a mapping to the 
two-dimensional antiferromagnetic Ising model on the trian- 
gular lattice, which was exactly solved by Wannier.i^E 

The underlying ensemble is special also in the sense that 
one does not only know the free energy and the entropy, but 
also the two-point correlation functions, at least asymptoti- 
cally. Since this is the autocorrelation of the system, when 
placing point scatterers of unit mass on each vertex point, the 
diffraction measure for almost all realisations of the lozenge 
random tiling (with edge length 1, say) is of mixed type, and 
has the form y = (y) pp + (y) ac - The pure point part is 17 

(r) PP = ^ E ((-i)V + (-i)* 2 P 2 +P 3 ) 2 W 

(23) 

where F* is the dual lattice of the triangular lattice, spanned 
by Vj = (l, — J-) and v 2 — (0, -^), and (k l ,k 2 ) * s a shorthand 
for the wave vector k l v 1 +k 2 v 2 6 r*. The pure point part 
reflects the underlying lattice structure, 5 while the absolutely 
continuous one is the fingerprint of the structural disorder. It 
is effectively repulsive in nature, as expected, which manifests 
itsel f 17 i 26 in the property that the diffuse intensity is 'repelled' 
by the Bragg peaks. 

The diffraction of the example from Figure |30]is shown in 
Figure [3TJ The pattern is lattice periodic. The pure point part 
(big spots) follows from the exact formula in (|23l , while the 
continuous part (small spots) was calculated numerically by 
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Hoffe 17 ' 64 via FFT techniques. 

The corresponding situation for the randomised Ammann- 
Beenker or octagonal tiling looks similar at first sight, and 
leads (via simpleton flip thermalisation, see Figure l32l to 
patches of the form shown in Figure [33] However, the pos- 
sible vertex positions are no longer restricted to a lattice, but 
only to the module Z[§ 8 ] with ^ g a primitive 8th root of unity. 
This module is the corresponding set of cyclotomic integers 
and a dense point set in the plane, as explained earlier. As a 
result, apart from the trivial Bragg peak at 0, the diffraction 
measure will be continuous, with singular and absolutely con- 
tinuous components. The reason behind this is the logarith- 
mically diverging fluctuation of the embedding surface from 
the deterministic surface of the model set relative. 62 Due to 
the larger positional freedom of the vertex points, this fluctu- 
ation is strong enough to destroy the coherence that is needed 
for non-trivial Bragg peaks, but not strong enough to avoid 
singular continuous contributions. 64 

Unfortunately, this is one of the claims that have not yet 
been proved, though there can be hardly any doubt about its 
correctness. A numerical calculation 16 of the diffraction of 
the finite patch shown in Figure [33] leads to the pattern of 
Figure [34] a similar result was obtained by Hoffe 64 via FFT. 




Fig. 32: A simpleton flip used in the thermalisation of the Ammann- 
Beenker (or octagonal) tiling. 



A comparison with the diffraction of the perfect Ammann- 
Beenker tiling in Figure [16] still reveals a lot of similarities, 
despite the approximative nature of the calculation. In partic- 
ular, one can clearly map the strong peaks of the perfect case 
to positions of the random tiling diffraction, and also various 
ring-type structures are clearly common to both images. In 
view of these similarities, it is not clear to what extent kine- 
matic diffraction of a finite patch can distinguish the perfect 
from a random tiling. 

The simpleton flip of Figure [32] provides a standard ap- 
proach for the preparation of random tiling samples. It works 
well also for other tilings with rhombic prototiles, where one 
might have different types of simpletons to consider (for in- 
stance, there are two such configurations in the rhombic Pen- 
rose tiling). One usually starts from a periodic approximant 
(to minimise boundary effects) to a perfect tiling, which is not 
difficult to construct, and runs the simpleton flip thermalisa- 
tion until correlations have decayed. In such ensembles, the 
process can be shown to be topologically transitive, so that the 
entire ensemble compatible with these boundary conditions is 
accessible.^ 4 - Note, however, that there are other important 
ensembles, such as the random square triangle tilings, where 
no such local flip exists. Here, one needs alternative methods, 
such as the well-studied 'zipper' moveJ- 02 that temporarily in- 
troduces some new (auxiliary) tiles that enable a randomisa- 
tion path, until the created tiles annihilate themselves again 
and leave a modified square triangle tiling behind. 

Finally, the case of random tilings in 3-space is clearly 
of great interest. A natural candidate from the very begin- 
ning^ 4 ^ 2 - has been the randomised version of the primitive 
icosahedral tiling, which is built from the two rhombohedra of 
Figure [18] While there are 24 complete vertex configurations 
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Fig. 34: Numerical approximation to the diffraction image of the 
random tiling of Figure [33] 

in the Kramer-Neri projection tiling, counted up to icosahedral 
isometries, there are 5450 possible ones in a typical random 
tiling. 30 So, it is clear that this version locally shows a much 
higher degree of disorder. However, unlike the previous exam- 
ples, the fluctuations away from the embedding hypersurface 
seem to be bounded, 44 i 62 which implies a diffraction of mixed 
type, this time with a pure point and an absolutely continu- 
ous component — despite the statistical icosahedral symme- 
try, which is non-crystallographic; a numerical confirmation 
was obtained by Monte-Carlo simulation techniques.^!! 

If one employs a statistical variant of the projection method, 
the fluctuations mentioned above lead to a distribution in in- 
ternal space that can be described by a density function. The 
latter will resemble a Gaussian profile, 44 ' 62 which makes the 
pure point part of the diffraction explicitly accessible via an 
appropriate extension of the model set theorem to this case. 1 1 1 
This gives diffraction formulas of PSF type where the sums on 
both sides run over dense point sets. A further generalisation 
was recently formulated for measures by Lenz and Richard. 90 

6 Outlook 

The discovery of quasicrystals 118 in 1982 had a profound im- 
pact on various disciplines, including mathematics and, in par- 
ticular, to harmonic analysis and mathematical diffraction the- 
ory. The approach described above emerged from the investi- 
gation of aperiodically ordered systems, and offers a method 
that can be applied to a wide range of structures. 

After 30 years of quasicrystal research, the diffraction 
of mathematical quasicrystals that are described by cut and 
project sets (regular model sets) is well understood. Such 



structures are pure point diffractive, much as conventional 
crystals, except that the Bragg peaks are supported on a point 
set that is dense in space. For many standard examples, the 
corresponding diffraction amplitudes can be calculated explic- 
itly, for instance in terms of Fourier transforms of the corre- 
sponding window(s). 

The situation changes quickly if one leaves the realm of 
model sets. As discussed, Meyer sets still inherit some of 
the structure, in the sense that their diffraction measure con- 
tains non-trivial pure point components. For substitution (or 
inflation) based structures, examples with all spectral types are 
known. In this article, we met examples of all three types — 
the Fibonacci chain (which is a pure point diffractive model 
set), the Thue-Morse chain (which has singular continuous 
diffraction) and the Rudin-Shapiro chain (with absolutely con- 
tinuous spectrum). In fact, it is easy to come up with a substi- 
tution system that has a mixed spectrum comprising all three 
spectral types. 

Quasicrystals are expected to contain some inherent (or 
structural) disorder, and it is therefore desirable to under- 
stand the effect of disorder on diffraction, and, vice versa, 
the conclusions on disorder that one can draw from examin- 
ing diffraction patterns, in particular with regard to continuous 
diffraction. This is far from being well understood, but the ex- 
amples discussed above provide a glimpse at the general situ- 
ation. As the Bemoullisation example shows, diffraction can- 
not always detect the nature of 'order', for instance whether 
the latter is of deterministic or entropic origin. Conversely, 
diffuse diffraction does not always need to be a sign of random 
disorder. At present, we only have a very limited knowledge 
of how large the homometry classes can be. In the pure point 
diffractive case, a recent approach by Lenz and Moody 88 ' 89 
provides one possibility for an abstract parametrisation. Un- 
fortunately, this approach does not seem to be extendable to 
cover continuous diffraction components. The investigation 
of further examples with different types or degrees of order 
will hopefully shed more light on this matter. 

One does not have to go far to find examples of important, 
yet still not completely understood systems. A prominent one 
is the Conway-Radin pinwheel tiling This tiling is based 
on a single triangular prototile (of edge lengths 1, 2 and V5), 
with an inflation rule of linear inflation multiplier y/5, so each 
re-scaled triangle (which is planar) is dissected into five con- 
gruent copies. Figure [35] shows a photograph of a patch of 
the tiling, which has been used as a theme for Melbourne's 
Federation Square development. Because the inflation con- 
tains a rotation that is incommensurate with n, a new direc- 
tion is introduced in each inflation step. Consequently, each 
infinite pinwheel tiling contains triangles in infinitely many 
distinct orientations, and the corresponding tiling space even 
has complete circular symmetry. 1Q0 i 106 i 107 The diffraction pat- 
terns shows striking similarity to a powder diffraction from a 
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Fig. 35: Detail of a facade at Melbourne's Federation Square featur- 
ing a pinwheel tiling. Photography © U. Grimm. 

square-lattice based structure. 7 While there is strong evidence 
for sharp rings in the diffraction pattern (which are singularly 
continuous in the plane), mimicking the case of the rotation- 
averaged square-lattice structure, the presence of further rings 
or absolutely continuous components is still unclear. 

More generally, one needs a unified setting for the diffrac- 
tion of systems with mixed spectra. An interesting suggestion 
was made by Gouere 53 on the basis of the intensity measure 
of the Palm measure of a point process. This provides an al- 
ternative way to define the autocorrelation of the system. It 
is possible to include cases such as crystallographic systems 
or model sets into this scheme, 88,89 and it was recently also 
shown 6 how to use this approach in a systematic way for sys- 
tems with various kinds of disorder. Since the theory of point 
processes is a highly developed branch 36 > 37 of modern proba- 
bility theory, the use of these methods looks rather promising. 
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